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~^ , Abstract 

In this paper, we establish the existence of a unique "regular" weak 
solution to the Large Eddy Simulation (LES) models of turbulence with 
critical regularization. We first consider the critical LES for the Navier- 
r~^ ] Stokes equations and we show that its solution converges to a solution 

of the Navier-Stokes equations as the averaging radii converge to zero. 
Then we extend the study to the critical LES for Magnetohydrodynamics 
equations. 
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1 Introduction 
> 

^b ' Let us consider the Navier-Stokes equations in a three dimensional torus T3, 

^ . divv = 0, (1.1) 

v,t + div{v(E)v) - iyAv + Vp = f, (L2) 

t^ ■ 

(~>j • subject to v{x,0) — Vo{x) and periodic boundary conditions. Here, v is the 

fluid velocity field, p is the pressure, / is the external body forces, i^ stands for 
the viscosity . 

Equations (|l.ip - (ll.2l) are known to be the idealized physical model to compute 
Newtonian fluid flows. They are also known to be unstable in numerical sim- 
. , ulations when the Reynolds number is high, thus when the flow is turbulent, 

j^ ■ Therefore, numerical turbulent models are needed for real simulations of turbu- 

lent flows. In many practical applications, knowing the mean characteristics of 
the flow by averaging techniques is sufficient. However, averaging the nonlin- 
ear term in NSE leads to the well-known closure problem. To be more precise, 
if denotes the filtered/averaged velocity field then the Reynolds averaged NSE 
(RANS) 

v.t + div(tJ (g) U) — vAv + Vp + div TZ{v, v) — f, (1.3) 



where TZ{v,v) — v iSi v — v (gi v is the Reynolds stress tensor, is not closed 
because we cannot write it in terms of 17 alone. The main essence of turbulence 
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modeling is to derive simplified, reliable and computationally realizable closure 
models. In |17| and |18| Lay ton and Lewandowski suggested an approximation 
of the Reynolds stress tensor, given by 



TZ{v,v) = V ®v — V <Siv. (1.4) 

This is equivalent form to the approximation 



div{v (>^ v) !v div{v 1^ v). (1-5) 

Hence, Layton and Lewandowski studied the following Large Scale Model con- 
sidered as a Large Eddy Simulation (LES) model: 

divw^O, (1.6) 

w^t + div{w (E)w) - lyAw + 'Vq=J, (1.7) 

considered in (0,r) x T3 and subject to w{x,0) = Wo{x) = vq and periodic 
boundary conditions with mean value equal to zero. Where they denoted {w, q) 
the approximation of (tJ, p) . 

The averaging operator chosen in (|1.7p is a differential filter, [TT], [T^], [S], [T7] . 
[9], [7], that commutes with differentiation under periodic boundary conditions 
and is defined as follows. Let a > 0, given a periodic function (p G £^(¥3), 
define its average Tp to be the unique solution of 

-a'^Alp + lp^ ip, (1.8) 

with periodic conditions, and fields with mean value equal to zero. 

The main goal in using such a model is to filter eddies of scale less than 
the numerical grid size a in numerical simulations. For a general overview of 
LES models, the readers are refered to Berselh et al. [5] and references cited 
therein. Notice that the Layton-Lewandowski model ()1.6|) - (|1.7p differs from the 
one introduced by Bardina et al. [4] where the following approximation of the 
Reynolds stress tensor is used: 



TZ{v, v) = V (El V — V 1^ V. (1-9) 

In [T7] and (TB] Layton and Lewandowski have proved that (ll.6p - (jl.7p have 
a unique regular solution . They have aslo shown that there exists a sequence 
aj which converges to zero and such that the sequence {Wa.,qa ) converges to 
a distributional solution {v,p) of the Navier-Stokes equations . 
We remark that many of these results established in the above cited papers 
have been extended to the following three dimensional magnetohydrodynamic 
equations (MHD): 

dtv - vxAv + div(t) ®v) - div(S ® ;B) + Vp 
dtB - V2AB + div(t; ®B)- div(S ® v) 

/ Bdx^ / t;da; = 0, div6 = divS 

6(0) = Z?o, 'y(O) = -yo, 



0, 


(1.10) 


0, 


(1.11) 


0, 


(1.12) 


^'0, 


(1.13) 
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here, v is the fluid velocity field, p is the fluid pressure, B is the magnetic field, 
and Vq and Bq are the corresponding initial data. The interested readers are 
referred to [HI [H] and references cited therein. 

This paper has two main correlated objects. The first one is to study the 
Large Eddy Simulation for the Navier-Stokes equations (LES for NSE) with a 
general filter — ^: 

divi(; = 0, (1.14) 



w^t + div(u; (g) It) ) — vAw + Vg = /, (1-15) 

c,2«(_A)V+^' = ¥', (1-16) 

where the nonlocal operator (— A)^ is defined through the Fourier transform 

(_AJV(fe) = |fepV(fc)- (1-17) 

Our task is to show that for > ^ (see Theorem 12. ip . we get global in time 
existence of a unique weak solution {w,q) to eqs. (|1.14p - (|1.16p such that: 
(it), q) are spatially periodic with period L, 

/ w{t,x)dx = and / q{t,x)dx = foriG[0,r), (1.18) 

JT3 JT3 

and 

w{0,x) =Wo{x) =v^ iuTg. (1.19) 

We note that the value 6* = i is optimal and of course the general a family 
considered here recover the case 9—1 studied in [TB]. The LES for NSE with 
9 = 1 can be also addressed as the zeroth order Approximate Deconvolution 
Model referring to the family of models in [T] . The value 6* = i is consistent 
with the critical regularization value needed to get existence and uniqueness to 
the simplified Bardina model studied in [3] . We note also that fractionnal order 
Laplace operator has been used in another a models of turbulence in [211 121 US] • 
The second object of this paper is to study the (LES) model for Magneto- 
hydrohynamics equation (LES for MHD) with a general filter — ^. Hence, we 
consider the following LES for MHD problem 



dtw ~ viAw + div{w (8) w ) - div{W ^W ) + Vg = 0, (1.20) 



dtW - V2AW + Ahj{w ®W )- dw{W ®w ) = 0, (1.21) 

/ Wdx= wdx = Q, div to = div VF = 0, (1.22) 

W{Q) = Wo.w{Q) = wo, (1.23) 

where the boundary conditions are taken to be periodic, and we take as before 
the same filter — ^ and {w^W^q) is the approximation of {v,B,p) solution of 
the MHD Equations. 
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The case when 6* = 1 is studied in [TS] where the authors gived a mathematical 
description of the problem, performed the numerical analysis of the model and 
verified their physical fidelity. 

In this paper, we show that for > ^ (see Theorem 13. ip . we get global in time 
existence of a unique weak solution {w,W,q) to eqs. (|1.20p - (|1.23p . Let us 
mention that the idea to consider the LES for MHD with critical regularization 
is a new feature for the present work. The Approximate Deconvolution Model 
for Navier-Stokes equations with > | is studied in ^ and the Approximate 
Deconvolution Model for Magnetohydrodynamics equations with 9 = 1 is stud- 
ied in [13] . As mentioned in [5] the value > | is not optimal in order to prove 
the existence and the uniqueness of the solutions in the deconvolution case. 
"The exponent "3/4" looks like a "critical exponent". We conjecture that we 
can get an existence and uniqueness result for lower exponents, but concerning 
the convergence towards the mean Navier-Stokes equations, we think that it is 
the best exponent, but this question remains an open one. " 
Notice however that unlike the LES case the value 9 = ^ is not sufficent to get 
the uniqueness in the deconvolution case. Based on this work, we will study 
in a forthcoming paper the Approximate Deconvolution Model for both Navier- 
Stokes equations and Magnetohydrodynamics equations with critical regular- 
izations. 

Finally, one may ask questions about the relation between the regularization 
parameter 9 and the model consistency errors. These questions are adressed in 
[TB] where 9=1. Therefore, the issue is to find the relation between the model 
consistency errors and the regularization parameter 9. 

This paper is organized as follows. In section 2 we prove the global existence 
and uniqueness of the solution for the LES for NSE with critical regularization. 
We also prove that the solution {w,q) of the LES for NSE converges in some 
sense to a solution of the Navier-Stokes equations when a goes to zero. Section 
3 treats the questions of global existence, uniqueness and convegence for the 
LES for MHD with critical regularization. 

2 The Critical LES for NSE 

Before formulating the main results of this paper, we fix notation of function 
spaces that we shall employ. 

We denote by ^^(Ta) and W^PiTs), s > -I, 1 < p < oo the usual Lebesgue 
and Sobolev spaces over T3, and the Bochner spaces C(0, T; X), LP{0, T; X) are 
defined in the standard way. In addition we introduce 

W^:'^ = Iwe VK'^P(T3)3; f ii; = 0; div ii; = in T3 j . 

We present our main results, restricting ourselves to the critical case 9 = ^, and 
for simplicity we drop some indices of 9 so sometimes we will write "^" instead 
of "^^", expecting that no confusion will occur. 



$LaTeX: 201 1/7/8$ 



2.1 Existence and uniqueness results for the LES for NSE 
Theorem 2.1. Assume that 9 = ^ 

i 2 

free function and Wq G W^^^ . The: 
solution to (|1.14p - (|1.19p such that 



Theorem 2.1. Assume that ~ ^. Let f G L'^{0,T;W^i- ) be a divergence 

- 2 

free function and Wq G W^Jiv ■ Then there exist {w,q) a unique "regular" weak 



weCiO,T;wt')nL\0,T;Wlt^''), (2.1) 

w^teL^iO,T;W-i'^), (2.2) 

gGL2(0,T;l¥^'2(T3)). (2.3) 



fulfill 



/ {w^t,<p) -{w(E)W,Vcp) + iy{Vw,Vcp) dt= {f,v)dt 
Jo Jo 

forallcpeL^O,T;wt^), 



(2.4) 



Moreover, 

w{0) = Wq. (2.5) 



Remark 2.1. The notion of "regular weak solution" is introduced in J^. Here, 
we use the name "regular" for the weak solution since the weak solution is unique 
and the velocity part of the solution w does not develop a finite time singularity. 

Remark 2.2. Once existence and uniqueness in the large of a weak solution to 
the model (|1.14p - (|1.19p with critical regularization is known. Further theoretical 
properties of the model with critical and subcritical regularizations can then be 
developed. These are currently under study by the author and will be presented 
in a subsequent report. 

Proof of Theorem 12. li The proof of Theorem 12.11 foUows the classical 
scheme. We start by constructing approximated solutions [v^ ,p^) via Galerkin 
method. Then we seek for a priori estimates that are uniform with respect to 
A''. Next, we passe to the limit in the equations after having used compactness 
properties. Finaly we show that the solution we constructed is unique thanks 
to Gronwall's lemma. 

Step l(Galerkin approximation). Consider a sequence {'^'^IJ^i consisting of 
L^-orthonormal and Ty^'^-orthogonal eigenvectors of the Stokes problem sub- 
jected to the space periodic conditions. We note that this sequence forms a 
hilbertian basis of L^ . 
We set 

JV JV 

«;^(t, x)^Y. ^r{tW{x), and g^(i, x) ^ Y. <lk (fV'''' ■ (2-6) 

r=l |fc| = l 
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We look for {w^ {t,x),q^ (t,x)) that arc determined through the system of 
equations 



(2.7) 
and 



Aq^ = -dwdw(n^{w^(g,w^)). (2.8) 



Where the projector H^ assign to any Fourier series \^ Qk^^^^ its N- 

fcGZ3\{0} 

dimensional part, i.e. 2. 9k^^ ^• 

fcGZ3\{0},|fc|<Ar 

Moreover we require that w^ satisfies the following initial condition 



N 



w^{0,.)^w^^J2''oV''{^), (2-9) 

and 

w^ ^Wq Strongly in M^S'2(T3)3 when iV -^ oo. (2.10) 

The classical Caratheodory theory [23] then implies the short-time existence 
of solutions to (|2.7l) - (|2.8p . Next we derive estimate on c^ that is uniform w.r.t. 
A'^. These estimates then imply that the solution of (|2.7l) - (|2.8p constructed on 
a short time interval [0,r^[ exists for all t G [0,r]. 

Step 2 (A priori estimates) Multilplying the rth equation in (|2.7p with 
a^^|fc|^^c^(i) + c^ (t), summing over r = 1, 2, ..., TV, integrating over time from 
to i and using the following identities 

-Aw^,w^ + aH^A)iw^) = \\w^\\l, + \\w^\\l^,^^^, (2.12) 

{7,w^ + Ji-A)iw'')^{f,w''), (2.13) 



and 



",V(u;^ + a*(-A)su;^)) = (ii;^ ® w^, Vii; 



divti) 



N liu 



(2.14) 



2 

leads to the a priori estimates 

^ J\f,w^) ds + ^[\\wo\\l + \\wo\\l 



(2.15) 
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Using the duality norm comined with Young inequahty we conelude from eqs. 
(|2?T5l) that 



sup H^||^+ sup \\w''\\l^ + ,, f (\\w''\\l, + \\w''\\l^^^ ds<C 
tG[0,T«[ te[0,T"[ '5' Jo ^ '^' ^ 

(2-16) 
that immediately implies that the existence time is independent of N and it is 
possible to take T = T^ . 
We deduce from p.l6p that 

w^ e L°°(0, T; W^l;^) n ^^(O, T; W^+^^^iTa)^). (2.17) 

From ((2Tf| and ([33]) it follows that 

rn;^eL2(0,T;VK^^2). (2.18) 



Consequently from the elliptic theory eqs (|2.8I) implies that 

i-T 



f Wp^'Wi.dtKK. 
Jo ''' 



(2.19) 



From eqs. dHZ]), (PTZ)) . ^^JE\i and ([^1^ we also obtain that 

i-T 







k^llia 2^^< -^- (2-20) 



Step 3 (Limit iV ^ oo) It follows from the estimates (PTf|) - (P^ and the 
Aubin-Lions compactness lemma (see |22| for example) that there are a not 
relabeled subsequence of {w^ , q^) and a couple {w, q) such that 

weakly* in L°°(0,r; VF^^^), (2.21) 

weakly in L^{0, T; W^+^'^), (2.22) 

wWw weakly in i2(0,r;W^S'2)^ (2 23) 

weakly in i2(o,r;VF"S'^), (2.24) 

weakly in i2(o,r;VF5'2(T3)), (2.25) 

strongly in L^{0, T; W'^^T^f), s < 1 + ^ (2.26) 

6 

wWw strongly in i2(0,T; VK''^2(T3)^),r< i (2.27) 

The above established convergences are clearly sufficient for taking the limit 
in (|2.7p and for concluding that the velocity part w satisfy (|2.4p . Moreover, 
from (|2.22p and (I2.24p one we can deduce by a classical argument ( see in [2]) 
that 

weC{0,T;W^''^). (2.28) 



7 







w^ 


-^* w 






w^ 


-^ w 


w^ 


(8) 


w^ 


^ w ^w 






w^, 


-^ Wj 






9^ 


^ q 






w"" 


— s> w 


w^ 


« 


w^ 


— > to (g) to 
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Furthermore, from the strong continuty of w with respect to the time with value 
in W^'^ we deduce that w{0) = Wq. 

Let us mention also that w+a^ [—A) ew is a. possible test in the weak formlation 
(|2.4p . Thus w verifies for all t e [0,T] the follwing equality 



\wit)\\l + \\wit)\\l,) +2:. 1'^ {\\w\\l, + \\w\\l 
\f,w)ds+(\\wo\\l + \\wo\\l 



(2.29) 



Step 5 (Uniqueness) Since the pressure part of the solution is uniquely 
determined by the velocity part it remain to show the uniqueness to the velocity. 

Next, we will show the continuous dependence of the solutions on the initial 
data and in particular the uniqueness. 

Let {wi,pi) and {w2,P2) any two solutions of (I1.14p - p.5p on the interval [0,T], 
with initial values Wi{0) and 11)2(0). Let us denote by 6w = Wi — W\. We 
subtract the equation for W\ from the equation for w^ and test it with Sw. We 
get using successively the relation p.5|) . the fact that the averaging operator 
commutes with differentiation under periodic boundary conditions, the norm 
duality, Young inequality and Sobolev embedding theorem: 

6 ' 6 ' 

< {w2 ®W2~'W\® wi,'V{Sw + a^{—A)^Sw)) 

< {W2 ^W2 — Wi ^ WijVSw) 



< -\\6w(g)Wi\\'^_i „ 



(2.30) 



4 

\X^..\\'^ ll».. II 

1+^ 



< -\\6w\\l J\w,''^ 



Using Gronwall's inequality we conclude the continuous dependence of the solu- 
tions on the inital data in the L°°{[0, T], VFS'^) norm. In particular, if Swq — 
then Sw — and the solutions are unique for all t G [0,T]. Since T > is 
arbitrary this solution may be uniquely extended for all time. 
This finish the proof of Theorem 12. II 



2.2 Limit consistensy for the critical LES for NSE 

In this section, we take the limit a — > in order to show the following result: 

Theorem 2.2. Let {wa, Qa) be the solution of (|1.14p - (|1.16p for a fixed a. There 
is a subsequence aj such that (Wa-,qa) ^ (f^TP) as j ^ 00 where {v,p) G 
L°°([0,r];L2(T3)3)nL2([o,r];W^]iV^)xii([0,T];ii(T3)) is a weak solution of 
the Navier-Stokes equations with periodic boundary conditions and zero mean 
value constraint. 
The sequence Wa ■ converges strongly to v in the space i^([0, T]; LPlT^)^) for all 
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2 < P < -^' ^^'^ weakly in L'^{0, T; L3r-4 (Ta)^) for all r > 2, while the sequence 
Qocj converges strongly to p in the space L^([0, T]; ^^(Ta)) for all -^ < p < ^, 
and weakly in the space L^ (0, T; L^^-'^ (T3)) for all r >2. 

Before proving Theorem 12.21 we first record the following three Lemmas. 

Lemma 2.1. Let 9 e R+,0 < /S < 29,8 eR and assume that ip e W^l^. Then 
Tp G W^div ' ■5'^c/i that 

II^IU+/3.2<4ll^lU,2, (2.31) 

and 

\m\s,2 < M\sa- (2.32) 

Proof, see in [2] 

Lemma 2.2. Assume Wa belongs to the energy space of solutions of the Navier- 
Stokes equations, then 

/ \\wa^ WaW's^^ dt < 00 for any -< r < 00. (2.33) 

Jo 2r .^ 6 

Proof. We have by interpolation that 

«;„eL'-(0,r;L3^(T3)3) (2.34) 

for any r > 2, thus we deduce by using Holder inequality that 

lOa^tCa eL5(o,r;L3l^(T3)3><3). (2.35) 

From Sobolev embedding we deduce that 

Wc,(g)Wc, ei5(o,T;l¥T^'2(T3)3x3)^ (2.36) 

for any r > |. 



Lemma 2.3. Assume Wa belongs to the energy space of solutions of the Navier- 

4 
3 



Stokes equations, then for all p > 1 and q > -k such that 



- + I- < 1, (2.37) 

p Sq 



have 



pT 3 + -3 

/ \\Wa®Wa-Wa®Wa\\ldt<Ca''''i' "' . (2.38) 

3p — 6 !y 

Proof. We take r = 2q, from the Sobolev injection W ^p ' (T3) ^-> L'^{T3), 
it is sufficent to show that 

/ \\wa(E)Wa -Wa<E)Wa\\%-e Jt< Ca'^^^p "'\ (2.39) 

Jo 2p ^^ 
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From the relation between Wa fX" Wa and Wa (S) w^ we have 

\\Wa ®Wa-Wa ®Wa\\'%^^ < a^'' \\w^rWw^\\lj^-6 „g „ (2-40) 



Lemma |2 . 1 1 imphes that 

Z 3r+2p-3pr t Z 

(2.41) 
Recall that 

/■"^ r 8 

/ ||i(;q, ® i(;q,|| i_3,, dt < oo for any - < r < oo. (2.42) 

This yields the desired result for any p > 1, g > | such that ^ + ^ < 1- 

Proof of Theorem 12.21 The proof of Theorem 12.21 follows the lines of the 
proof of the Theorem 4 in ^8] . The only difference is the strong convergence 
of the pressure term qa to the pressure term p of the Navier-Stokes equations. 
We will use Layton-Lewandowski [IF, as a reference and only point out the 
differences between their proof of convergence to a weak solution of the Navier- 
Stokes equations and the proof of convergence in our study. First, we need to 
find estimates that are independent from a. Using the fact that Wa belong to 
the energy space: L°°{[0,T]; L^iT^f) D L^{[0,T];W^l^) and from the Aubin- 
Lions compactness Lemma (the same arguments as in section 2.1) we can find 
a subsequence {wa-,qa) and {v,p) such that when aj tends to zero we have: 

w^^ -* V weakly* in L°°i[0,T]; L^T^f), (2.43) 

w^^--v weakly in i2([o,r];W^i^2(T3)3), (2.44) 

Wc,^ -^ V weakly in i''(0, T; L^^ (Tg)^) for all r > 2, (2.45) 

Wa^ (g) Wa^ -^ V (g) V weakly in L^{0,T;L^^ {Tsf^), for aU r > 2, (2.46) 

ti;„^ ^ V strongly in ^^([0, T]; ^^(Ta)^) for all 2 <p < y, 

4 5 

Wa^ ® Wa, -> f ® v Strongly in LP{[0, T]; ^^(Tg)^''^) for aU - < p < - 

o o 



(2.47) 
5 

3' 

(2.48) 



Having (12.38^ and (|2.48p at hand we deduce that 

T3)3x3) foraU - <p< 

(2.49) 



mr^ ^v^v strongly in LP(0, T; LPCIsf^) for aU - < p < -. 

3 O 3 



Then, from ([235]) and (|2^ we deduce that 
Wa^ (gWa^ --v(g)v wcakly iu L^ (0, T; L^fei (Tg)^''^) for all r > 2. (2.50) 

10 
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Further, we have 



qit) ^ R{J2w'^^w^t)) (2.51) 

k,l 



where the linear map R is defined by 

R : L'iTsf ^^ L%T3) (2.52) 

{u''')k,i=i,2,3 ^ i-Ay'dkdiiu''') (2.53) 

By the theory of Riesz transforms, i? is a continuous map for any s e]l,oo[. 
Consequently, from (I2.50p we have 

/ Ika II '3,. dt < 00, for all r > 2. (2.54) 

Jo ' ^^^ 

From p.49p we deduce that for almost alH > 0, 

4 5 

Wa^ ^Wa^t) -^v^v{t) strougly iu LP(T3)^^^ for all - < p < - . (2.55) 

o o 

Using the dominate convergence theorem and the continuity of the operator R, 
we conclude that 

q^^ ^ p strongly in ^^([0, T]; ^^(¥3)) for aU ^ < p < | (2.56) 

Finally we deduce from (|2.56p and (I2.54p that 

<?a, ^P weakly in L^ {0,T;L^^{T3)), for ah r > 2. (2.57) 

These convergence results allow us to prove in the same way as in [TBI that 
{v,p) is a weak solution to the Navier-Stokes equations, so we will not repeat 
it. 

3 Application to the LES for magnetohydrody- 
namic equations (LES for MHD) 

In this section, we consider the critical LES regularization for magnetohydro- 
dynamic (LES for MHD) equations, given by 



dtw - vi Aw + div {wSi W) -div{W ^W) + \/q = 0, (3.1) 



dtW - V2.AW + div(t(; ®W^)- dm(yV ® lo^ ) = 0, (3.2) 

/ Wdx^ wdx^O, divi(; = divVF = 0, (3.3) 

W{0) = Wo, w{0) = Wo, (3.4) 



11 
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where the boundary conditions are taken to be periodic, and we take as before 
the same spacing average operator. 

a3(-A)6(^6 + (p6 = (y5, ip{t,x + Lsj) = ip{t,x). (3.5) 

Here, the unknowns are the averaging fluid velocity field w{t,x), the averag- 
ing fluid pressure q{t,x), and the averaging magnetic field W{t,x). Note that 
when a = 0, we formally retrieve the MHD equations. Existence and unique- 
ness results for MHD equations are established by G. Duvaut and J.L. Lions 
in [To]. These results are completed by M. Sermange and R. Temam in [20] . 
They showed that the classical properities of the Navier-Stokes equations can 
be extended to the MHD system. 

The aim in this section is to extend the results of existence uniqueness and 
convergence established above for the LES for NSE to the LES for MHD. We 
know, thanks to the work [15], that for 9 = 1 these results hold ture. Further, 
when = i , we proved in the above section the existence of a unique "regular" 
weak solution to the LES for NSE. Therefore, it is intersecting to find the critical 
value of regularization needed to establish global in time existence of a unique 
"regular" weak solution to LES for MHD. 

We divide this section into two subsections. One is devoted to prove the 
existence of a unique "regular" weak solution to the LES for MHD with = ^■ 
The second one is devoted to prove that this solution converges to a weak 
solution to the MHD equations when a tends to zero. 

3.1 Existence and uniqueness results for the LES for MHD 

First, we establish the global existence and uniqueness of solutions for the LES 
for MHD equations with = ^. 
We have the following theorem: 

1 - 2 

Theorem 3.1. Assume that 9 = ^. Assume Wo and Wo are both in W^^]^ . 
Then there exist (w,W,q) a unique "regular" weak solution to (|3.ip - p.4p such 
that 

w, W^C{Q,T-wl^)rM\Q,T-wlt^-^), (3.6) 

w^u W ^t ^ L\Q,T;W-i-^), (3.7) 

q^L^iQ.T-.w'^^-'iT^)). (3.8) 

fulfill 

I {w,t,(p} - (wWw, Vip) + (WWW, Vif) + vi{Vw, Vip) dt = 
Jo 

/ (W,f , tp) - {w(g)W, V(p) + {W (E)w, Wip) + V2{Vw, V(p) dt = 
Jo 

forallipeL^iO,T;wt^). 
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Moreover, 

u?(0) = it-o and W{Q) = Wn- (3.10) 

Proof, of Theorem 13.11 We only sketch the proof since is similar to 
the Navier-Stokes equations case. The proof is obtained by taking the inner 
product of pr^ with a^{-A)^w + w, (fr2T1) with ai{-A)iW + W and 
then adding them, the existence of a solution to the critical LES for MHD can 
be derived thanks to the Galerkin method. Notice that {w,W) satisfy the 
following estimates 

[\\wit)\\i + \\wmi^) + [wwmi + wwrni^) 

= \\\wo\\l + ||«;o|li J + (llWolll +iw^o|li J ■ 

(3.11) 
The averaging pressure q is reconstructed from w and W (as we work with pe- 
riodic boundary conditions) and its regularity results from the fact that w (^ w 



3x3\ 



It remains to prove the uniqueness. Let {wi, Wi, qi) and {w2, W2, 92), be two 
solutions, 6w = W2 — Wi,SW — W2 — Wi, Sq = q2 — qi- Then one has 



dtSw — viASw + div(u;2 (^ W2) — A\m{wi ® Wi) 



- dw{W2 ® W2) + dw{Wi ® Wi) + V5q = 0, 
dtSW - V2MW + div(«;2 ® W2) - div(toi ® Wi) 



- d:\NiyV2 ® W2) + div(VFi <S) wi) = 0, 

and Sw = 0, SW = at initial time. One can take a3[—A)eSw + 6w as test 
in the first equation of p.l2p and a^{—A)e6W + SW as test in the second 
equations of (|3.12p . Since Wi is divergence- free we have 

/ / wi(g)dw -.VSw ^ - / {wi-y)Sw -Sw^O, (3.13) 

Jo Jts Jo Jts 

Thus we obtain by using the fact that the averaging operator commutes with 
differentiation under periodic boundary conditions 

rT f 

(div(i(;2 ® W2) — div(t(;i ® if 1)) • ( a^ (— A) e (Jio + bw 

i-T /. 

(div(ti;2 ® 1^)2) — div(K;i ® Wi)) ■ Sw (3.14) 

^ f 

/ 5w W2 : VSw. 

Jts 

Similarly, because (if 1) is divergence- free we have 

/ I wi®SW : VSW = - / I {wi- V)SW -SW ^0, (3.15) 

Jo Jtr Jo JT:, 
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and thus we have the following identity 

i-T 



I I (div(u;2®W2)-div(u;i®VFi)) • (a^(-A)5(5W + (5W^) 

= f f {dW{w2(»W2)-dW{wi(»Wi))-SW (3.16) 

Jo JTa 



f 5wi»W2: V(5W. 

JT3 



Concerning the remaining terms we get by integrations by parts and by using 
the using the fact that the averaging operator commutes with differentiation 
under periodic boundary conditions 

f I {-diY{W2®W2)+diY{Wi®Wi)) ■ (a^{^/^)-^6w + 5w\ 

= 11 {~d\Y{W2®W2) + diY{yVi®Wi))-5w (3.17) 

Jo JTs 

f Wi(^SW : S/5w + f f SW®W2: VfJiw. 
Jtr Jo Jts 



and similarly 



/ / {- div {W 2 <E) W2) + div {W 1 (g) w 1)) ■ (a^-A)hw + SW'] 

«/ " Tl 3 

= f [ idiv{W2®W2)~div{Wigwi))-SW (3.18) 

Jq JT3 

= -/" [ {Wi-V)5w5W+ [ [ SW(x)W2:VSW. 
Jo JT3 Jo Jt3 

Therefore by adding (|3.14p - (|3.18p and using the fact that the averaging op- 
erator commutes with differentiation under periodic boundary conditions we 
obtain 

-^ / (\Sw\^+ai\vhwf)+^( [ \dW\^ + a^\\/hw 
2"^ JT3 ^ ^ 2ai \Jt3 

\''T3 / \>/ T3 

5w ®W2: V5w + Sw(g W2 : VSW 



Ta __ JT3 

SW(dW2: S76w - / SW ® W2 : VSW. 



(3.19) 
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By the norm duality 

1/ (5u;(8)t(;2 : V(5i(;| < ||(5ti;®t(;2ll_i 2l|V(5u;||i 2, (3.20) 

t [ Sw(^W2: V5W| < \\6w W2II 1 2l|V(5Ty|| 1 2, (3.21) 

/ 6W (E)W 2 ■■\75w\ < \\6W (E)W2\\_i oW'^SwWi ^, (3.22) 

I 5W®W2:V5W\<\\5W®W2\\_i2\W5W\\i2. (3.23) 
A3 



/T; 

By Young's inequality. 



— \\5W®W2f_i ^ <- 

V2 f ■ V2 



Hence, 



5u;||i 2 + \\5W\n 2 
n, 



1/ 5w®W2:V5w\<—\\5w®W2\\'ti2 + ^-\\^^w\\l2^ (3-24) 

JT3 ^1 6^ 4 6' 

1/ 5w®W2:V5W\<—\\5w®W2f_io+'^\\^5W\\l,, (3.25) 

7t3 ^^2 6'^ 4 6.^ 

I / (51^ (g) PF2 : V(5t(;| < —\\5W ® W2f-_i , + ^||V(5u;||i „ (3.26) 

Jt3 vi 6^^ 4 6.-^ 

1/ 5W®W2:V5W\<—\\5W®W2f_i^+'^\\y5W\\l^. (3.27) 

JT3 ^2 6- A 8' 

By Holder inequality combined with Sobolev injection 



— \\5w®W2f_i^ < — ll'^u'fi Jk2||?, 1 2 (3.28) 

±\\Swr^W2\\\, < -\\Sw\\lJW2\\l^, , (3.29) 

-l!<5Ty®W2f_i2 < -¥w\\ljW2\\l+. ^ (3.30) 



— ||(5W®-u;2f_i,2 < — ll'5W^lli,2ll«'2||?+i,2 (3.31) 



+ min(j/i,z^2) (||'5u?||?,2 + ll'5W^lll2) +'^^ min(i^i,i^2) (||'5u?||^^i 2 + ||(5W^||?+i,^ 

< —7^ 7 (\\Sw\\l , + \\SW\\l ^ (\\W2\\1^, 2 + i|W-2||?+i 2) 

mm(i^i,i^2) ^ "5 15'^ \ '-^6'^ -^+6'^/ 

(3.32) 
We conclude that du = 6B = thanks to Gronwall's Lemma. 
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3.2 Limit consistensy for the critical LES for MHD 

Next, we will deduce that the LES for MHD with critical regularization gives 
rise to a weak solution to the MHD equations. 

Theorem 3.2. Let {wa,Wa,qa) be the solution of (|1.14p - (|1.16p for a fixed 
a. There is a subsequence Uj such that (Wa ,Wa ,qa ) -^ {v,B^p) as j —> oo 
where iv,W,p) G [L-([0, T]; L2(T3)3) n L^qq, T]; V]^)]^ x lI ([0, T]; ii (T3)) 
is a weak solution of the Navier-Stokes equations with periodic boundary condi- 
tions and zero mean value constraint. 
The sequence Wa converges strongly to v in the space LP([0,T]; LPiT^)'^) for 

all2<p< ^, and weakly m L''(0, T; L^^ (Ta)^) for all r > 2. 
The sequence W a^ converges strongly to B in the space LP{[0, T]; LPiT^)^) for all 
2 < p < -y , and weakly in L''(0, T; L^'--^ (Ta)^) for all r > 2, while the sequence 
qa- converges strongly to p in the space L^([0, T]; LP(T3)) for all -^ < P < §, 
and weakly in the space L^ (0, T; L'^-^-^ (Ta)) for all r > 2. 

Proof of Theorem 13.21 As in the proof of Theoreni l2.2l we can show that 
for all I < p < I we have 

strongly in LP{0,T; LPiTsf"^), (3.33) 

strongly in 1^(0, T; LPiTsf"^), (3.34) 

strongly in LP{0, T; LPiTsf''^), (3.35) 

Wa(S)Wa^B^B strongly in LP(0,r;iP(T3)^''^). (3.36) 

The above LpLp convergences combined with the fact that w^ and W^ be- 
long to the energy space of the solutions of the Navier-Stokes equations and 
the Aubin-Lions compactness Lemma allow us to take the limit a ^' in order 
to deduce that {wa, Wa,qa) converge to {v, B,p) a weak solution to the MHD 
equations. The rest can be done in exactly way as in |18] , so we omit the details. 

Acknowledgement: The author thanks professeur R. Lewandowski for 
interesting discussion about this paper. 
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